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Abstract 

Using a recently developed semi-analytical method (Self-Consistent Local RPA or SC-LRPA) 
we study the stability of the ferromagnetic phase in diluted magnetic systems where the exchange 
coupling between magnetic impurities are of RKKY form. A short discussion of the relevance of 
these calculations with respect to the ferromagnetism observed in diluted ferromagnetic materials 
is provided. Then, within a two step approach, we study ferromagnetism in Zni^xCr^Te. In 
the first step of our study, we calculate the magnetic couplings between Mn impurities within the 
LDA. In the second step, we diagonalize the resulting effective Heisenberg Hamiltonian using the 
SC-LRPA. We also compare, when available, our calculations with Monte Carlo simulations and 
experimental measurements. 
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Diluted Magnetic Semiconductors materials have attracted considerable attention be- 
cause of their technological potential in the domain of spintronics. These materials are 
particularly promising since a small concentration of magnetic impurities can already give 
rise to a relatively high Curie temperature . The key issue to obtain quantitative and 
reliable values of the Curie temperature is twofold. First, one has to be able to provide 
correct values for the exchange couplings between magnetic impurities. Secondly, one must 
treat properly the resulting effective Heisenberg Hamiltonian which describes the interac- 
tion between randomly distributed magnetic impurities on the lattice. By properly, we mean 
provide an accurate treatment of (i) thermal fluctuations and (ii) effects of disorder in the 
position of the impurities. Different approaches have been developed for the determination 
of the exchange integrals in Ill-V materials. There are essentially two different kind of 
methods. First, simple phenomenological approaches 0, |^ in which the band structure of 
the host material is described in a realistic manner by a Kohn-Luttinger Hamiltonian for 
the valence band states and a coupling Jpd (or Vpd) between the itinerant carrier and the 
localized magnetic impurities. In these approaches, the coupling between localized magnetic 
impurity and the itinerant carrier is only treated perturbatively. However, within a simple 
single band model in which the local Jpd coupling was treating non perturbatively within 
CPA, it has been shown that the perturbative treatment which leads to Jij oc J^^ and thus to 
Tq oc is only restricted to very weak couplingy|. A second class of approaches are based 
on first principle calculations. In these approaches, the coupling between d states of the 
magnetic impurities and host states is treated non-perturbatively, and this leads to realistic 
values of the exchange couplings. The main, and great, advantage of these approaches is 
that the couplings do not depend on adjustable parameters. 

This manuscript is organized as follows. The first part is devoted to the study of the 
effects of thermal fluctuations and disorder on the stability of the ferromagnetic region as the 
density of itinerant carriers in diluted magnetic systems, with assumed RKKY-like exchange 
couplings. In view of our results, we discuss the relevance of such a phenomenological model 
to simulate the exchange couplings in actual diluted ferromagnetic semiconductors. In the 
second part, using realistic first principle calculations, we calculate the exchange couplings 
within the Tight-Binding LMTO approach in Cr doped ZnTe semiconductors. We then 
calculate the Curie temperature as a function of the magnetic impurity concentration for 
this material. 
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Throughout this paper, we present calculations of the Curie temperature obtained us- 
inga recently developed Self- Consistent Local Random Phase Approximation (SC-LRPA) 



j. The Curie temperatures calculated within SC-LRPA appear to be in very good a gree - 
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ment with Monte Carlo simulations where the same exchange couplings have been used 
However the SC-LRPA has several important advantages with respect to Monte Carlo calcu- 
lations: it allows (i) to derive a direct semi-analytical expression for the Curie temperature 
(no extrapolation is required), (ii) the spins can be either quantum or classical, (iii) the 
calculations are very fast and (iv) correlation in the disorder can be included without any 
complication |2| . Within SC-LRPA the self-consistency is usually achieved after approxi- 
mately 10-15 iterations, i.e. it requires the same number of diagonalizations of a Ni^p x A^imp 
complex matrix. For comparison, for a given configuration of disorder, our calculations of 
Tc are at least 3 orders of magnitude faster than Monte Carlo calculations. Of course speed 
is not the only point; in the low temperature phase the SC-LRPA includes quantum effects 
that cannot be simulated classically. We now summarize the main steps of the SC-LRPA 
treatment. The disordered Heisenberg Hamiltonian which we consider is 

The sum runs only over pairs of sites occupied by magnetic impurities. The spins are 
quantum with value S. In the case of classical spins, we will properly perform the limit 
5 — s> oo at the end of the calculations. We now define the following retarded Green's 
function, 

/•-I- DO 

oo 



G,,{uj) = / G,,{t)e'^'dt (2) 



where G,,{t) = -id{t){[St{t)-,Sjm)- 

After standard Tyablicov decoupling of the equation of motion of Gij{uj) we obtain, 

{u - ht")G,,{u) = 2{St)5,, - {St) E (3) 

I 

where the local effective field is, 

K'f = Y.Mst) (4) 
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For a given temperature and fixed disorder configuration, the local magnetization (Sf) 
has to be determined self-consistently at each impurity site. In order to close the set of 
equation we use Callen expression, which relate the local Green's function at site i to the 
local magnetization at this site (iflj. 



1 + - $2 



The local effective magnon occupation number reads, 

(h = -1 /'+°° ImGujuj) 
' 27i{Sf) J-oo exp{u/kT) - 1 ^ ^ 

The previous set of equations allows us now to determine, at each temperature, the local 
magnetization at each impurity site, and the dynamical properties such as the dynamical 
structure factor. 

To determine the value of the Curie temperature we take the limit {S?) ^ in the 
previous set of equations. We obtain, 

kBTc = \siS + l)-^j:y. (7) 

"J imp i 

where, 

(8) 

The local spectral function Aii{E) = —^Im{ ^"^^^ ). In the previous equation the factor 
\i = liranp^Tc^^ (^^ is the magnetization averaged over all impurity sites). It is interesting 
to note that the Curie temperature can be re-expressed in term of the eigenf unctions and 
eigenvalues of the following effective Hamiltonian, where the matrix elements are. 



I 



This lead to. 



F.^E^ (10) 



E 
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Thus the Curie temperature becomes, 



= ^s{s + 1) E(E ^^1^)-^ (11) 

'J-' 'imp j ct Q 

This expression imphes that the nature (extended/locahzed) of the eigenstate will have 
important effects on the magnetic excitation spectrum and on the Curie temperature. A 
detailed discussion of these aspects will be provided elsewhere jll|. 

We now turn to the commonly used model of dilute magnetic systems where we place 
spins randomly, and independently, on a host lattice, and assume an interaction between 
each pair of spins of RKKY form with some density of carriers mediating the exchange. 
The carrier density is an important parameter, as it determines the scale of oscillation of 
the exchange. In Fig.l, we have plotted the Curie temperature for a fixed concentration 
of magnetic impurities as a function of the carrier density, normalized with respect to the 
concentration of impurities. We assume that the exchange couplings are of RKKY form and 
we including, the effects of disorder in an averaged way, i.e. by including an exponential 
damping of the couplings. Thus the interaction is taken to be 

r (sin(2A;i;'r) - 2A;i;'r cos(2A;i;'r)) , . 

- '^oexp(--^) (12) 

with Fermi vector kp determined by the density kp = (Svr^nc)^ where Uc = is the density 
of carriers for impurity concentration x. 

The curves of Curie temperature versus carrier concentration are shown for three different 
values of damping. Note that we do not consider in this work the effect of a short range 

n 

super-exchange contribution; this is discussed in detail in Ref. |5j. For each case we show, for 
comparison (dotted lines), the Curie temperature estimated from mean- field Virtual Crystal 
Approximation (MF-VCA). This approximation has often been used, but, unlike ours, re- 
places the random lattice of impurities by a translationally invariant effective medium. We 
observe that as the carrier concentration increases from zero, our calculated Curie Temper- 
atures have very different behaviour from the mean-field theory. Most dramatically, rather 
than increasing monotonically with carrier density, the Curie temperature saturates and 
vanishes at a density of carriers much smaller than the density of impurity moments. This 
is significant in that in doped Ga(Mn)As, for example, it is clearly seen from simultaneous 
transport and magnetic measurements that the highest Curie temperatures are obtained 




FIG. 1: Curie temperature as a function of the carrier density for a fixed concentration of magnetic 
impurities (3% ) and exchange couphngs of RKKY form (see equation 12). The exponential cut-off 
in interactions is defined by the screening parameter . 

as the samples are annealed, giving carrier densities per moment approaching unity. Thus 
we conclude that simple RKKY interactions are not appropriate to describe the ferromag- 
netism in doped III-V semiconductors. The reasons for this had been obscured by the use 
of mean-field theories. The reason that the SC-LRPA shows the instability at higher carrier 
concentration is that includes the effect of the oscillating tail of the exchange constants, i.e. 
the frustration in the magnetic interactions. As damping becomes more short range, i.e. Ro 
becomes smaller, the frustrating tail of interactions is reduced, as seen in Figure 1, but this 
is not enough to suppress the instability. 

We learn from this that in order to explain observed ferromagnetism in diluted magnetic 
semiconductors we need to model the magnetic interactions more accurately. An approach 
that has proved successful, both in doped III-V Q and II- VI materials, as we shall see below, 
is to use ab initio methods to estimate the magnetic couplings at all distances including 
effects of the (disordered) band structure, correlations, screening treated within, for example 
Local Density Approximation (LDA). As the magnetic couplings can be calculated by a 
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FIG. 2: Calculated Heisenberg exchange couplings (in Kelvin) as a function of the distance (in 
units of the lattice constant) between magnetic impurities for different concentration of Cr x = 3% 
, 5 %, 10 % and 20 % . 

magnetic force theorem, i.e. in principle from the energies of the ground state perturbed by 
weak magnetic fields, they should be accurate. With the resultant Hamiltonian, which is 
specific to both the host material and the doping concentration (since disorder and doping 
effects are taken into account in the band structure), we will apply our method of SC-LRPA 
to estimate the Curie temperature. Within the class of II- VI semiconductors we turn now 
to the experimentally promising candidate for room temperature ZnTe doped with Cr. 

First we show the exchange integrals between Cr moments deduced from our ab initio 
calculations. In Fig. 2, we have plotted the magnetic exchange integrals between Cr 
impurities in ZnTe as a function of the distance, for different concentration. First, we 
clearly see that the frustration, which, we argued, could strongly suppress ferromagnetism 
in an RKKY model, is absent here. This is related to the fact that the Chromium 
dopants are close to the configuration Cr^"*", i.e. do not, to a first approximation, introduce 
carriers into the II- VI host. In fact the band structure shows "half-metallicity" : a finite 
density of states at the Fermi energy for the majority carriers and a gap for the minority 
carriers [l^. The couplings, while oscillating in magnitude, remain essentially ferromagnetic 
and exponentially damped. The ferromagnetic nature at short distances is attributed to 
the strong hybridization of the Cr 3d levels with the Te 5 p states [14]. The damping has 
different origins, the effects of disorder (treated within CPA), the half-metallicity and the 
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effects of the d levels. As the concentration of Cr increases, the couplings decrease at short 
distances keeping the same position of local minima. 
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FIG. 3: Calculated Curie temperature within Self-Consistent Local RPA and Mean Field VCA for 
Zni^xCvxTe as a function of x. 



In Fig. 3 we show the variation of the Curie temperature calculated within SC-LRPA as a 
function of the magnetic impurity concentration. A systematic average over at least 100 con- 
figurations of disorder has been taken for each concentration. For comparison, we also show 
the Mean-Field VCA results. While qualitatively the differences are not as dramatic as in 
Fig.l in that both mean- field and our curve are monotonically increasing, quantitatively the 
predictions are very different. Most importantly, we believe our theory should be a reliable 
guide to the highest Curie temperatures that should be obtained for each concentration Q]., 
In particular we predict that room temperature should require concentrations around 25 %; 
the mean field VCA gives much lower values for this minimum concentration, about 2%. 
Experimentally Saito et al^^ have found = 300 3 10 for thin films with x ~ 0.2 which 
is close to our prediction(T(7(x = 0.2) = 265). Monte Carlo methods jl^ for Zn(Cr)Te us- 
ing the same exchange couplings were made for a few concentrations and are close to our 
curve. For example the Monte Carlo result is (Tc(x = 0.2) 300). We would conclude that 
the samples measured are close to the optimal Curie temperature for the concentration of 
uncorrelated impurities. 

In conclusion, we have shown that the region of stability of ferromagnetism in diluted 
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systems where the exchange couphngs are RKKY-hke is very narrow even when there is large 
damping. These couphngs are very different from those obtained in ferromagnetic diluted 
semiconductors by means of first-principle calculations, and thus appear to be inappropriate 
to explain the wide region (in carrier concentration) of stability usually observed in these 
materials. We have calculated the variation of the Curie temperature in Zni_xCr,j.Te as a 
function of Cr density and found that room temperature can be achieved for densities above 
X ~ 0.25. This is important for this particular material and also illustrates how, by coupling 
first principle calculations of the exchange couplings with SC-LRPA, we may make useful 
predictions for a range of important new materials. 



[1] H. Ohno, Science 281,951 (1998). 

[2] K. W. Edmonds et al, Phys. Rev. Lett. 92,037201 (2004), K. W. Edmonds K.Y. Wang,R.P. 
Campion, B.L. Gallagher, C.T. Foxon, Appl. Phys. Lett. 81,4991 (2002). Additional values of 
Tc were provided by Edmonds ct al ( private communication). 

[3] T. Dietl, A. Haury and Y. Merle d'Aubigne, Phys. Rev. B 55, R3347 (1997) 

[4] M. Abolfath, T. Jungwirth, J. Brum, and A. H. MacDonald, Phys. Rev. B 63, 054418 (2001) 

[5] R. Bouzerar, G. Bouzerar and T. Ziman submitted to Phys Rev B. 

[6] G. Bouzerar, J. Kurdnovsky and P. Bruno, Phys. Rev. B 68, 205311 (2003) 

[7] G. Bouzerar, T. Ziman and J. Kudrnovsky, Eur. Phys. Letters, 69, 812-818 (2005). 

[8] G. Bouzerar, T. Ziman and J. Kudrnovsky, Phys. Rev. B 72, 125207 (2005). 

[9] G. Bouzerar, T. Ziman and J. Kudrnovsky, Appl. Phys. Lett. 85 4941 (2004). 
[10] H. B. Callen, Phys. Rev. 130,890 (1963). 
[11] G. Bouzerar et al., in preparation. 

[12] L. Bergqvist, O. Eriksson, J. Kudrnovsky, P.A. Korzhavyi, and I. Turek, Phys. Rev. Lett. 93, 
137202 (2004). 

[13] H. Saito,V. Zayets, S. Yamagata and K. Ando , Phys. Rev. Letters 90, 207202 (2003) 
[14] Q. Wang, Q. Sun,P. Jena and Y. Kawazoe, J. Appl.Phys. 97 043904 (2005). 



9 



